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Abstract 

The structure function, S(k), for the s = 1, Haldane gap antiferromagnetic 

chain, is measured accurately using the recent density matrix renormalization 

group method, with chain-length 100. Excellent agreement with the nonlinear 

a model prediction is obtained, both at k ~ it where a single magnon process 

dominates and at k ~ where a two magnon process dominates. We repeat 

our calculation with crystal field anisotropy chosen to model NENP, obtaining 

good agreement with both field theory predictions and recent experiments. 

Correlation lengths, gaps and velocities are determined for both polarizations. 
75.10.-b, 75.10.Jm, 75.40.Mg 



1 



Typeset using REVIgX 



It was first argued by Haldane Jl| that integer-spin s, antiferromagnetic chains exhibit a 
gap to a triplet magnon above a singlet ground state. This is by now well confirmed both 
experimentally |2]|| and numerically ■ Much insight into these systems can be obtained 
from the nonlinear a (NLa) model, the long-distance field theory limit. Since this model 
only becomes an exact representation for large s, it is a priori unclear how accurate it is 
for s = 1, the case of greatest experimental interest. Furthermore, the NLcr model itself 
is not very tractable and is therefore often approximated by a weakly interacting or free 
boson (Landau-Ginsburg) model |J. The NLcr model gives information about two limiting 
wave-vectors: k — > it and k — > 0. The excitations near it are predicted to be a completely 
stable single-magnon, plus a three (and higher) magnon continuum. Near the excitations 
are a two (and higher) magnon continuum. The single-magnon contribution to the structure 
function has a simple square root Lorentzian (SRL) form near k ~ 7r which only depends 
on the magnon interactions via a renormalization of the overall scale. On the other hand, 
near k « the form of S(k, us) depends quite strongly on the interactions, but the overall 
scale is fixed by a sum rule since the total spin operator must obey the spin commutation 
relations. One non-trivial exact result is known about the NLcr model which is of relevance 
to possible experiments on spin chains: the explicit form of the two-magnon contribution to 
S(k,u) at small k [jT0|l . The NLcr model does not predict how S(k,u) crosses over from one 
magnon to two magnon behavior as k is swept across the Brillouin zone. Anisotropy may 
be included in the NLcr model, but this destroys the integrability and we must then use the 
free boson approximation to the two-magnon contribution. 



Recently, detailed inelastic neutron scattering experiments |TT[ were carried out on 
Ni(C 2 H 8 N 2 )2N0 2 (C104) (NENP), measuring S(k,u) for k > .3vr. At all k, the peak width 



appeared to be resolution limited. With a ratio of inter- to intra chain coupling estimated fT2 
to be 4 x 10 -4 , NENP is a highly one dimensional compound. It is, however, not isotropic. 
The largest of the anisotropies is the single ion anisotropy, D. Neglecting other smaller 
anisotropies NENP is well described by the Hamiltonian 
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H = J^{S ? ,-S m + J D(^) 2 }, (1) 



with D/ J = 0.18 and J = 3.75meV [13|. Calculations of S(k, u) using exact diagonalization 
for L < 18 < 16 @) or Monte Carlo for L < 32, (g have found some evidence for 

the two-particle nature of the small k excitations from the presence of spectral weight above 
the lowest-lying state. 

Here, we present accurate results on the equal time structure function S(k), using the 
density matrix renormalization group (DMRG) method recently proposed by White ||, with 
chain-length, L = 100. Details of our calculation will be given elsewhere [jl(J. Our results 
for S(k) are similar to those published previously for D/J = 0.0 (L = 18) |§, D/J = 0.18 
(L = 16) |15| except for a significant reduction of finite-size effects near k = tt and the 
presence of data at smaller k. This improved precision is crucial for much of our analysis. 
We compare our results on the pure Heisenberg model to the NLcr model predictions. The 
SRL form is obeyed accurately for k > .8n, with a large renormalization of the amplitude. 
The two-magnon NLcr model result, including the overall amplitude, is remarkably accurate 
for k < Arc. Using previous Monte Carlo results |J for u> e (k), the minimum energy excitation 
at wave-vector k, we show that the single mode approximation (SMA) to S(k) is very good 
for k m it, indicating that the three (and higher) magnon contribution is small. However, 
the SMA fails at small k suggesting the two-particle nature of the excitations there. We 
also obtain results for a single-ion anisotropy term with D/J = 0.18 to model NENP. We 
determined precise values for the singlet and doublet gaps, A.n/J = 0.6565(5), A±/J = 
0.2998(1) (in good agreement with Hl5]), correlation lengths, £|| = 3.69(5), £j_ = 8.35(7) 



and velocities, v\\/J = 2.38(1), v±/ J = 2.53(1), and the ground state energy per spin 
e /J = —1.2856861. These obey the relativistic relationship, £ a = hv a /A a to within 2% for 
both polarizations. Again the SRL form is very accurate near k ~ 7r and the free boson 
model is fairly accurate at small k. The SMA is very accurate at large enough k but fails 
badly at small k. Good agreement is obtained with the experiment apart from an overall 
scale factor of about 1.25 which was to be expected since the experimental total intensity 
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exceeded the exact sum rule: (1/L) Y, a kS aa (k) = s(s + 1) by 30% (± 30%). The behavior 
of S(k) suggests that slightly higher resolution and lower k neutron scattering experiments 
may reveal the two-particle nature of the excitations at small k. 

We implement the DMRG using density matrices of size 243 x 243 keeping 81 eigenvectors 
at each iteration. Our results are obtained using a finite lattice method on a chain of length 
100. For a discussion of the numerical procedure we refer the reader to the discussion by 
White ||. The DMRG yields a higher accuracy if open boundary conditions are used, and 
all of our results are therefore obtained for open ended chains. If subject to open boundary 
conditions the s = 1 Heisenberg antiferromagnetic chain has a singlet ground state with 
positive parity, + , and an exponentially low lying triplet, 1~ [[L7],|18|]. In the thermodynamic 
limit the ground state is thus four fold degenerate. We have performed our calculations in 
the 1~ state since this is computationally more convenient. We have repeated some of the 
calculations for the + state and for a smaller system of 60 site. In none of these cases 
were the results seen to differ aside from boundary effects which were seen to diminish with 
increasing chain length. 

Let us first discuss our results for the isotropic chain where D = 0.0 in Eq. (|I]). The 
gap A/J = 0.4107 @J3JT§ and the velocity v/J = 2.49(1) |H| is known. Using the DMRG 
method we calculated the bulk correlation function < S% S? >, from which the correlation 
length £ = 6.03(1) can be determined [|l^] in agreement with previous results 0. Since our 
calculations are performed in the 1~ state < Sf > has a non zero expectation value giving 
rise to a disconnected part in the correlation function. The bulk correlation function we 
consider is therefore < S% Sf > — < Sj >< Sf >. We find that end effects are so small 
that we may assume this correlation function to be equal to the correlation function for a 
periodic chain. The static structure factor, S(k), can now be calculated by a simple Fourier 
transform S aa (k) = J2r ex P('>'kr)S aa (r), where S aa (r) is the above equal time correlation 
function in real space. Our results for S = S xx = S yy = S zz are shown in Fig. [I] as the open 
squares. The results obey the sum rule (1/L) J2 a k S aa (k) = s(s + 1) to within numerical 
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precision up to small boundary effects. 

As already mentioned the isotropic s — 1 Heisenberg antiferromagnetic chain can be 
described by the NLcx model UJ, with Lagrangian density 

where a is the lattice spacing and describes the sublattice magnetization. We have ap- 
proximately Si = (— l) l s<p + al, with 1 = l/(vg)4> x (d(j>/dt). Since we have set fi = 1, the 
velocity, t> , has dimension energy times length. The massive triplet of fields 4> is restricted to 
have unit magnitude, <fi 2 = 1. The lowest energy excitation corresponds to a single magnon 



with k = 7r and a relativistic dispersion -E(fc) = y (A) 2 + {v) 2 {k, — n) 2 . Near k = ir our 
numerical results for have a maximum corresponding to this stable single magnon, 

exhibiting the relativistic SRL form: 

QV 1 (3) 



2 y/A 2 + i; 2 (A; - tt) 2 

Excellent agreement between this form and the numerical results is seen for k/n > 0.8 (see 
Fig. p]). Since v and A are known from independent calculations this allows us to determine 
the coupling constant, g ~ 1.28. This result can be compared to the usual large s result of 
g = 2/s = 2 and tentatively to the 1/s expansion result |19| g ~ 1.44 for the bare coupling. 

Near k = S(k) approaches zero as expected since the ground state is a singlet. In 
this region S(k) is expected to be dominated by a two (and higher) magnon continuum 
corresponding to the excitation of two magnons with momenta ±7r. Exact results for the 
two magnon contribution to S(k, u) within the NLcr model frame work have been obtained 



in Ref. JlC ] : 



v k 2 Juj 2 - (vk) 2 - 4A 2 

Here uo 2 — {vk) 2 = 4A 2 cosh 2 (#/2) and |G(^)| 2 is given by the expression 

\C(ftW 2 - ^i+M. / tanh(g/2) \ 2 



This result can be numerically integrated over to to yield the two magnon contribution to 
the static structure factor. The result is shown as the solid line in Fig. [I] and, as can be 
seen, it is remarkably accurate out to k/n < 0.4. 

If the nonlinear constraint, (p 2 = 1, on the field in the NLcx model is lifted, a simpler 
model of almost free bosons can be arrived at || with Lagrangian density 

By allowing the 3 velocities and gaps to be phenomenologically adjusted, anisotropy can be 
included. If the A</> 4 term is neglected a mean field theory of free bosons is obtained. In the 
isotropic case, within this free boson model, the two magnon contribution to S(k) can be 
shown [[H|] to be equal to Eq. (^) with |G(#)| 2 = 1. The free boson estimate is shown as the 
long dashed line in Fig. [1]. While in good qualitative agreement with the numerical results 
near k = the discrepancy at larger k, between this result and the exact NLa prediction, 
demonstrate the importance of interaction effects that are neglected in the free boson theory. 

Within the SMA the dynamical structure factor S(k, u) is approximated by a 5-function 
at the single mode frequency c^sma, S(k,u) = So(k)6[u — cjsma(^)]- Since the first moment 
of the dynamical structure factor obeys a sum rule pi this implies |TB| that 



to, 



sma(*0 = -J(F X + F y )(l - cos(k))/S zz (k), (7) 



with F a =< S^S°" +1 >. Denoting the lower edge of the spectrum by o> e , one can show |2l|JT6 



that the inequality cjsma(^) > oJ e (k) is obeyed. Furthermore it can be shown [0 that if 
S(k,u) has a 5-function peak at the frequency ou§(k), which is below the lower edge of the 
continuos part of the spectrum, u c (k), then u>sMA{k) > 0Js{k). This inequality breaks down 
if u)g(k) > u c (k). Using the numerical results F a = e /3 = —0.4671613 for the isotropic 
chain we can now calculate cusmaW from S(k). Our results are shown in Fig. along with 
the quantum Monte Carlo (QMC) results of Takahashi |J for the lower edge, u e , of the 
spectrum. Also shown is the relativistic dispersion relation (solid line), which should be 
very accurate near k = it, and the bottom of the two-magnon continuum (long dashed line), 
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2y A 2 + v 2 {k/2) 2 , valid near k = 0. As can be seen in Fig. ^]both predictions coincides with 
the QMC results. While highly accurate near k = it the SMA is clearly seen to break down 
near k = 0. From our results for the structure factor we expect the two magnon scattering 
to become important around k ~ 0.47T which seems to be consistent with the divergence 
of wsma from the QMC results. If compared to the relativistic dispersion, which should 
represent the single magnon dispersion relation ujk near k = it very accurately, it is seen 
that also here there must be a small multi- magnon contribution since c^sma is above Uk- 

We now turn to a discussion of our results for the anisotropic chain with D/J = 0.18. 
Using the DMRG we have calculated < S§ S^ > and < S£ Sf > in the 1~ state. In this 
case the disconnected part, < S% >< >, is so small that we can neglect it. Fourier 
transforming we obtain S" = S zz and S 1 - = S xx = S vy , respectively. Our results are shown 
in Fig. ^. We estimate end effects to give an error of 0.05% in S^lk = it), somewhat larger 
at k — 0. End effects are judged to be negligible for S 1 ". Also shown in Fig. are the INS 



results of Ma et al. [it which are directly comparable to our results. The open triangles are 
points where to within experimental accuracy 5*" and S 1 - were identical. The full squares 
and circles are data points where S" and S"- 1 , respectively, could be resolved experimentally. 
At k = 7r the experimental results are about 20-30% higher than our numerical results. 
Taking into account a 25% [^] uncertainty in the overall (multiplicative) normalization of 
the experimental results we find a good agreement between the numerical and experimental 
results. 

Allowing for different velocities, gaps and coupling constants for the two modes we expect 
S aa (k) to be well described by the SRL form Eq. (0) near k = it. With the gaps and velocities 
mentioned above and with g\\ = 1.44, g± = 1.17 this SRL form is shown as the dashed lines 
in Fig. [3| Very good agreement when k > 0.8tt between the SRL form and the numerical 
results is evident for both S" and S- 1 . 

In the region near fc = 0we see that since rotational symmetry around the x and y axis 
is broken, S 1 - can take on a non zero value at k = 0. This is clearly seen in the numerical 
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results for S . We can extend the results of Ref. JT0| for the free boson prediction to allow 



for different velocities. For the two magnon contribution to S aa (k) near k = we obtain (16 

2 f dk'dk" u b k ,y c u c k „v b 



with cyclic permutation of the indices a, b, c. Here u k = yA^ + (v a k) 2 . These results are 
shown as the solid lines in Fig. ^. As was the case for the isotropic chain the free boson 
prediction seems to be qualitatively correct near k = 0. The crossings of S^(k) and S ± (k) 
around k = 0.857T and 0.17T (see Fig., are reproduced correctly by the free boson model. 
In the range 0.1% < k < 0.857T S*" is the larger of the two structure factors. 

We have repeated the SMA calculation for the anisotropic case finding qualitatively 
similar behavior as for the isotropic chain. Again the SMA fails for small k while giving 
relatively precise results near k = 7T, indicating a small but non zero contribution from the 
multi magnon continuum at k = n. 

We would like to thank Steven White for sendings us Refs. |6]|?j prior to publication 
as well as for useful discussions. We are also indebted to S. Ma, C. Broholm, D. H. Reich, 
B. J. Sternlieb, and R. W. Erwin for permission to present their data. ESS thanks Dan Reich 
for several helpful discussion and for suggesting looking at the single mode approximation 
for this problem. This research was supported in part by NSERC of Canada. 
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FIGURES 

FIG. 1. The structure factor, 5, vs. k/ir for a 100 site isotropic chain. The numerical results 
are shown as open squares. 

FIG. 2. Dispersion of excitations vs. k/n for a 100 site isotropic chain. The points are QMC 
data from Ref. || for L = 32. The short dashed line is cjsmaW obtained from S(k) for L = 100. 

FIG. 3. S 1 " and S 1 - vs. k/ir for a 100 site anisotropic chain. The numerical results are shown 



as open squares and circles, respectively. The data points are INS data from Ref. [11|. 
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